Band-structure of graphene above the vacuum level 
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We study the 2D band-structure of the stand-alone graphene in the context of the bound states 
immersed in the 3D continuum. Depending on the symmetry of the in-plane k vector, the true 
bound states either survive in the continuum or not. By use of a nearly exactly solvable model we 
demonstrate that those bound states which do not survive acquire a finite life-time, in other words 
they turn into resonances. 



The low-lying energy bands of a stand-alone graphene 
are well understood. The linear dispersion of its tt (high- 
est occupied) and tt* (lowest empty) bands can be ob- 
tained analytically in the tight-binding approximation 
and it is known since as early as 1947 Presently, 
it is an easy exercise in the application of any of the ex- 
isting solid-state codes to calculate the tt as well as the 
lower cr valence bands from the first principles [ll-Q] . 

Going higher in energy, we eventually enter the contin- 
uous spectrum which corresponds to the infinite motion 
perpendicular to the layer. This is presented in Fig. (TJ 
where the band-structure of graphene obtained in the 
repeated-geometry calculation is plotted. Inside the con- 
tinuous spectrum one can discern the lines entering 
from below the continuum. The latter have a clear inter- 
pretation: They correspond to the states of the in-plane 
motion but with the energy above the continuum edge 
0. If, however, one looks deeper into the problem, an 
important issue arises: An electron placed in a state of 
the in-plane motion with a sufficiently high energy has, 
generally speaking, a non-zero probability to escape from 
the layer into vacuum, which would impart this state a 
finite life-time, i.e., turn it into a resonance. 

Since the discreet states inside the continuum, both 
true bound states and resonances, would play an im- 
portant role in optical absorption and inverse photoc- 
mission of graphene, to mention only two experimental 
techniques, they deserve a close study. The latter is the 
subject of the present work. It must be also noted that 
although our discussions are centered around graphene, 
the significance of the problem under investigation ex- 
tends beyond the latter, being pertinent to any stand- 
alone 2D-periodic system. 

We start by considering a trivial case: Let us have a 
quantum well Vz{z) in z direction with the flat poten- 
tial in xy plane. Then, if the well has bound states, and 
since the two perpendicular motions are independent, the 
wave-function is the product of a bound state in z direc- 




FIG. 1. (color online) The band-structure of graphene 
obtained in the repeated-geometry. The all-electron full- 
potential linearized augmented-plane wave code Elk [l^ was 
used for this calculation. The separation between the period- 
ically stacked layers is d = 400 bohr. 



tion and a plane wave in xy plane. As a result, there exist 
states which arc bound within the well while having an 
arbitrarily high energy above the vacuum level due to the 
motion in xy plane. 

If, however, we apply a potential which is periodic in 
xy plane, the variables in the Schrodinger equation do 
not separate any more, i.e., the two perpendicular mo- 
tions become coupled. To get a better insight on how 
this affects the high-lying energy bands of the in-plane 
motion, we first consider a simplified model: 

S-function quantum well with laterally periodic poten- 
tial ~ We are looking for a solution of the Schrodinger 
equation 
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with the model potential being a product of a periodic 
function in the xy plane and the (5-function quantum well 
in z direction 



y(z,r„) = ^VGe'°-"-a(z), 
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where G are the 2D reciprocal lattice vectors. We set 
Vo < to ensure the existence of a state bound to the 
z — Q plane. The solutions of Eq. ^ with the potential 
([2]) can be written explicitly as Bloch waves with respect 
to the motion in xy direction 
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where k is the in-plane wave-vector within the first Bril- 
louin zone, and ac are still unknown coefficients 
Importantly, in Eq. ([3|) we have retained the exponent 
with one sign only, which corresponds to considering the 
bound and, possibly, resonant states, while omitting the 
scattering states with respect to the z direction motion 
12| . The jump in the wave-function's z-derivative is ob- 
tained by the integration of Eq. ([T]) in z over the infinites- 
imal interval [0_,0+] 

V''(0+, r„ ) - V^'CO- , r II ) 2 ^ ^ce'^-^ii V(0, r , ), (4) 
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where the prime denotes the differentiation with respect 
to z. Together, Eqs. ^ and (g]) lead to 

^ ^(G - G')aG' = iV'^E - (G + k)2 aa- (5) 

G' 

The crucial point is the choice of the sign of the square 
roots in Eqs. ^ and ([5]). Denoting the generic square 
root by s, the rule must be 

if Re > then Re s > 0, else Im s > 0, (6) 

which choice ensures the correct asymptotic behavior of 
both bound states and resonances. 

The values of E which allow for non-trivial zero solu- 
tions of the homogeneous system of linear equations ([5]) 
determine the band-structure of our model system. How- 
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ever, in contrast to the original Eq. i 
tute a nonlinear eigenvalue problem 
E. This complication is the price paid for the separation 
of the bound and resonant states from the continuous 
spectrum. 

In Figs. [2]and[31 results of the numerical solution of the 
non-linear eigenvalue problem ([S]) which give the band- 
structure of the stand-alone layer are presented together 
with the results of calculations carried out for the same 
system in the repeated-geometry. The symmetry of the 
2D periodic potential was chosen that of a honeycomb 
lattice with grahene's lattice constant, and the values of 
the Fourier coefficients of the potential in atomic units 



FIG. 2. (color online) Band-structure of the model system 
obtained with the repeated-geometry calculation (blue lines) 
and with solving the eigenvalue problem ([5]) for a stand-alone 
plane (red points) . The calculation has been conducted along 
T - K - M - V line. 
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FIG. 3. (color online) The same as Fig. O but for the asym- 
metric k direction along 7bi + 13b2, where bi and b2 are the 
primitive reciprocal vectors. 



were: Vo = -1.093, Vb = 0.156 (-7/10 6 and 1/10 &, re- 
spectively, where h is the length of the primitive vector 
of the reciprocal lattice) for the first G star, and Vg = 
otherwise. In Fig. [2l the wave- vector varies along the 
F — X — M — F lines while in Fig. |3]a completely asym- 
metric direction of k variation is chosen. For the bands 
below the vacuum edge, both the repeated-geometry and 
the stand-alone calculations yield the identical results re- 
gardless of the symmetry of the wave- vector. In contrast, 
above the vacuum edge, whether or not a particular state 
localized near z = Q survives as a true bound state is de- 
termined by the symmetry of the corresponding k point. 
For the totally asymmetric case of Fig. [31 there are no 
such states. As can be seen from Fig. [21 along the sym- 
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TABLE I. Eigenenergies (in eV) for the model system ob- 
tained with the reduced size of V{G — G') matrix permitting 
the fuUy analytical solution of the non-linear eigenvalue prob- 
lem ©. 
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metric directions some of the bound state bands do sur- 
vive. Moreover, an isolated high-lying bound state exists 
in the F point. 

While the results in Figs. [5] and [3] elucidate the issue of 
the bound states survival above the vacuum edge, they do 
not answer the question what happens with those which 
do not survive, i.e., do the latter turn into resonances by 
acquiring finite life-times or they disappear at all. This 
is due to our numerical search for the eigenvalues of the 
nonlinear eigcnproblcm ([5]) having been restricted to the 
real axis of since no decisive numerical procedure ex- 
ists to cither find all complex- valued roots of this problem 
or to prove their absence. To shed the light on the latter 
issue, we have also been solving the eigenvalue problem 
([5]) analytically for the sizes of the V{G — G') matrix ad- 
mitting such a solution within a reasonable time. Using 
the Mathematica symbolic algebra software, we have been 
analytically evaluating the determinant A{E) of the sys- 
tem ([SJ, consecutively eliminating the square roots from 
the equation A.{E) — 0, which permitted its reduction 
to the polynomial equation. All roots of the polynomial 
equation (including the complex ones) have then been 
found with no loss of any of them guaranteed. Since spu- 
rious zeros have been largely introduced when reducing 
the equation to the polynomial, the roots were finally 
sorted retaining only those which satisfied the original 
equation A{E) = 0. This has been done for F, K, and 
M points with the matrices sizes of 19, 13, and 7, respec- 
tively. While all true bound states, both below and above 
the vacuum edge, were found to reproduce those previ- 
ously obtained numerically, in addition, complex eigen- 
values, which correspond to resonances, have been found. 
Results of this calculation are collected in Table HI 

Resonances in the excitation spectra of a quantum- 
mechanical system physically manifest themselves in clas- 
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FIG. 4. (color online) Energy dependence of the coefficient 
of reflection of an electron from our model system. The in- 
cidence is normal to the layer (k = 0). The vertical dashed 
lines indicate the position of the resonances at F point from 
Table HI 
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FIG. 5. (color online) z-dependence of the wave-functions of 
the states above the vacuum level in graphene obtained in ab 
initio repeated-geometry calculation. Plots a) and b) suggest 
bound states inside the continuum at F point and a point 
on the V — K line. In contrast, plot d) suggests a resonant 
character of the states on the T — M line. Plot c) exemplifies 
the wave-function of an arbitrary continuous spectrum state. 

tic scattering at real energies in the vicinity of the com- 
plex eigenvalues [l^l ■ In Fig. 0] we plot the coefficient 
of reflection of an electron incident normally onto our 
model system. The features in the reflectance spectrum 
clearly agree with the resonances' positions listed in the 
first column of Table H] (F point). The advantage of us- 
ing the reflection spectrum for studying the resonances is 
that this task does not require the solution of the eigen- 
value problem and it is performed at real values of the 
energy. 

From the model, we now return to the graphene prob- 
lem. The only information on the character of the states 
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above the vacuum level that can be derived from an ab 
initio calculation within the repeated-geometry is the lo- 
calization of the wave-functions around the layer. In 
Fig [SI we plot z-dependence of the modulus squared 
of the wave-functions of some representative states in F 
point, on the T - K, and on the T - M lines [ll]. For 
the chosen states in F point and on the T — K line, the 
wave-function behavior suggests a bound state with no 
detectable probability of an electron to be found in the 
intermediate region between the layers. On the contrary, 
for the points on the F — M line, the probability spills 
out of the layer and has a nonzero distribution over the 
whole of the inter-layer region. This behavior suggests 
that the state is not a bound one although the spill-out 
of the electron density is small. Our previous model anal- 
ysis suggests that these states become resonances having 
acquired final life-times. The third from the left inset in 
Fig. [5] shows the wave-function of an arbitrarily chosen 
state of (quasi-) continuous spectrum. 

In conclusion, we have found that the band-structure 
of a stand-alone 2D-periodical quantum-mechanical sys- 
tem above its vacuum level is strongly affected by the 
possibility for a particle to escape from the correspond- 
ing state into vacuum. The overall picture of the 2D 
bound states' either survival or disappearance has been 
drawn using a simple model of a periodic potential in 
xy plane complemented with the (5-function well in the 
perpendicular direction. Our analysis performed within 
this model suggests that the possibility for a particle to 
escape into vacuum leads to the bound states of low sym- 
metry turning into resonances, while the bound states of 
higher symmetry survive up to high energies above the 
vacuum level. 

Our model analysis has proven relevant to the ab ini- 
tio characterization of high-lying 2D energy-bands of the 
stand-alone graphene confirming that the true bound- 
states either survival or turning into resonances is deter- 
mined by the symmetry of the corresponding k point. 
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